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Abstract 

We investigate supergravity models in which the lightest supersymmetric particle (LSP) is a 
stable gravitino. We assume that the next-lightest supersymmetric particle (NLSP) freezes out 
with its thermal relic density before decaying to the gravitino at time t ~ 10 4 — 10 8 s. In contrast 
to studies that assume a hxed gravitino relic density, the thermal relic density assumption implies 
upper, not lower, bounds on superpartner masses, with important implications for particle col- 
liders. We consider slepton, sneutrino, and neutralino NLSPs, and determine what superpartner 
masses are viable in all of these cases, applying CMB and electromagnetic and hadronic BBN con- 
straints to the leading two- and three-body NLSP decays. Hadronic constraints have been neglected 
previously, but we find that they provide the most stringent constraints in much of the natural 
parameter space. We then discuss the collider phenomenology of supergravity with a gravitino 
LSP. We find that colliders may provide important insights to clarify BBN and the thermal history 
of the Universe below temperatures around 10 GeV and may even provide precise measurements 
of the gravitino's mass and couplings. 

PACS numbers: 04.65. +e, 12.60.Jv, 26.35. +c, 98.80.Es 
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I. INTRODUCTION 



Supersymmetric theories predict the existence of a spin 3/2 particle, the gravitino, the 
partner of the spin 2 graviton. The gravitino mass is 

F 

where F is the scale of supersymmetry breaking, and M* = {SitG n)~ 1 / 2 — 2.4 x 10 18 GeV is 
the reduced Planck mass. The masses of scalar superpartners are derived from terms such 
as 

Ay / d'O^p^ , (2) 

J 1V1 mcd 

where Ay are unknown constants, Z is a superfield whose auxiliary component develops the 
vacuum expectation value F, $j are standard model superfields, and M me( j is the mass scale 
of the interactions that mediate supersymmetry breaking. Similar terms give the spin 1/2 
superpartners mass. In supergravity, the mediating interactions are gravitational, and so 
M me d ~ M*, F ~ (10 10 GeV) 2 , and the gravitino and all standard model superpartners have 
mass ~ F/M* ~ M wca k, with the precise ordering determined by unknown constants, such 
as Ay. 

Most studies of supergravity have assumed, either explicitly or implicitly, that the lightest 
supersymmetric particle (LSP) is a standard model superpartner. This avoids potential 
complications resulting from the decay of standard model superpartners to a gravitino LSP, 
which naturally happens at time t ~ 10 4 — 10 8 s, well after Big Bang nucleosynthesis (BBN). 
However, the phenomenology and cosmology of gravitinos have also been considered in a 
number of studies El IH El H IH IE H, IM 1^ liol Ull fl2l] . (See also related studies of axino and 
quintessino dark matter [13J, [l4j, UJJ [l6( . ) 

More recently, it has been shown that the gravitino LSP possibility does not destroy the 
beautiful predictions of BBN even when the gravitino LSPs produced in late decays have 
relic density fi n = 0.23 and so are present in sufficient numbers to account for all of dark 
matter 17], |l8|, [l9| . In fact, bounds from the cosmic microwave background (CMB) and, in 
some corners of parameter space, entropy production and the diffuse photon spectrum may 
be even more severe than bounds from BBN (l7L 0] . Nevertheless, all of these bounds were 
shown to be respected for some regions of parameter space with weak scale superpartners. 
The possibility of superweakly-interacting massive particle (superWIMP) gravitino dark 
matter from NLSP decays thu s appears to be viable. The analogous scenario in extra 
dimensional theories 
related scenarios 3, 
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18, 2(|, as well as interesting astrophysical implications in this and 
22| have also been discussed. 
In this work, we take an approach that differs from the exploration of superWIMP grav- 
itino dark matter. Instead of assuming that gravitinos are the dark matter with Qq = 0.23, 
we assume that the NLSP reaches its thermal relic density ^nlsp De f° re decaying, and so 
Qq = (^g/mNLSp)^NLSP- Th & t i s ! we relax the constraint that gravitinos from NLSP de- 
cays account for all of dark matter. Rather we assume the simplest thermal history for the 
Universe and ask what regions of (^^,mNLSp) parameter space are allowed. The thermal 
relic density assumption has consequences that differ markedly from the fixed gravitino relic 
density assumption. To see this, assume that the gravitino and NLSP masses are both 
parametrized by a general superpartner mass scale msusY- As noted in Ref. |l9j |. if one 
assumes a fixed gravitino relic density, the NLSP number density scales as 1/msusY- Low 
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superpartner masses are therefore disfavored. In contrast, if one assumes a thermal relic den- 
sity for the NLSP, fift LSP oc (cxv}- 1 

m susY! where (av) is the thermally-averaged NLSP 
annihilation cross section. The NLSP number density then scales as mgusY, and so high 
superpartner masses are disfavored. This difference has obviously important implications 
for collider searches for new physics, and we discuss collider implications below. 

Even given the NLSP thermal relic density assumption, gravitinos may still be all of 
dark matter — for example, if the gravitino relic density from NLSP decays is too low, the 
remainder may be made up by gravitinos produced during reheating. However, the exis- 
tence of such alternative gravitino sources is either untestable, or testable only with strong 
assumptions about the early Universe. In contrast, the existence of a gravitino compo- 
nent from NLSP decays makes several robust predictions that are testable at cosmological 
observatories and collider experiments, and we concentrate on this gravitino source here. 
Before leaving the topic of reheating altogether, however, we note that the gravitino LSP 
scenario has an important virtue with respect to reheating. For stable weak-scale grav- 
itinos, the overclosure constraint is well-known to require a bound on reheat temperature 
of Tr < 10 10 GeV jnj. Recently, however, it has been shown that if the gravitino is not 
the LSP, hadronic BBN constraints greatly strengthen this bound [23|. For example, if the 
gravitino decays to the LSP + hadrons with branching fraction 10 -3 , the reheat tempera- 
ture must satisfy Tr ^ 10 6 GeV. This is uncomfortably low. The gravitino LSP scenario 
is therefore preferred if one requires a high reheat temperature, as might be desirable, for 
example, for scenarios of leptogenesis 

In the present analysis, in addition to the constraints on electromagnetic ener gy r elease 
considered previously, we include the recent results on hadronic BBN constraints 2^L The 



work of Ref! represents a significant update to previous hadronic analyses B0HII 
To include these results correctly, we must, of course, determine the leading contributions 
to hadronic energy. For slepton 1 and sneutrino NLSPs, the leading contribution is from 
three-body decays 

I -> IZG , vWG 

z> -> vZG , IWG . (3) 



The three-body decays have been studied in Ref. [19[- For a neutralino NLSP, the leading 
contribution to hadronic energy is from the two-body decays, such as 

X - ZG, hG , (4) 

followed by Z, h — > qq. These decays, and the hadronic constraints on them, were neglected 
in previous works. As we will see, however, they are the leading constraints in much of 
parameter space and they are especially important when the superpartner masses and their 
splittings are all of the order of the weak scale, the most natural possibility. 

The decays of Eqs. (|5J) and (JH) may be suppressed kinematically if mNLSP — m c < m z, mw 
or dynamically, as when the neutralino is photino-like. However, even in these cases, decays 
such as I — > IqqG and x ~~ * qqG are still possible at higher order. We have included 
estimates of these in our analysis. These decays are in some sense "model-independent"; 
even in the extreme case where the dominant decay is to invisible particles, at higher order 
there will be contributions to hadronic cascades from such decays. The hadronic bounds 



Throughout this work, "slepton" refers to a charged slepton. 
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are so constraining that these should be considered for any late decaying particle, whether 
a superpartner, an axino, a modulus, or other particle. 

After determining the regions of parameter space allowed by cosmology, we discuss the 
collider signals. The upper bounds on superpartner masses resulting from the thermal relic 
density assumption imply promising prospects for superpartners to be within reach of future 
collider experiments. In addition, we will see that the signals of supersymmetry in gravitino 
LSP scenarios may be completely different from the conventional supersymmetry signals. In 
particular, if sufficient NLSPs can be collected and monitored for decays, the NLSP lifetime 
may be measured, which may considerably sharpen our understanding of BBN and the 
thermal history of the Universe at temperatures of 10 GeV and below. Such studies may 
also provide the first direct measurements of the gravitino mass and the Planck scale from 
particle physics [29, 30]. 

The gravitino LSP possibility, assuming a thermal NLSP relic density, has been discussed 
recently in the context of minimal supergravity jjy]]. Our work is complementary in that we 
do not work in a specific model framework, but rather consider several NLSP candidates, 
as well as a wide range of gravitino and NLSP masses. Our work also differs in that we 
consider the hadronic constraints and the leading two- and three-body decays that contribute 
to hadronic energy. As noted above, we find that these are the leading constraints in the 
most natural regions of parameter space. 



II. LATE DECAYS 

We first discuss the decays of NLSPs for each of the various NLSP cases. NLSPs freeze 
out and are highly non-relativistic when they decay. We will be most interested in deriving 
the electromagnetic (EM) and hadronic energy releases 

6 = ei^NLSP , (5) 

where i = EM, had, because BBN constraints are, to a good approximation, constraints on 
£em and £had- Here B x is the branching fraction into EM/hadronic components, and e; is 
the EM/hadronic energy released in each NLSP decay. These are discussed in this section. 
y NLSP = n NLSP /n^ G is the NLSP number density just before NLSP decay, normalized to 
the background photon number density n® G = 2£(3)T 3 /7r 2 . 2 Given the assumptions of this 
work, Inlsp is determined by the thermal relic density for each NLSP; it is discussed in 
Sec. ED 

We will consider the cases of slepton, sneutrino, and neutralino NLSPs. As specific 
examples in each of these categories, we will focus on f R , u T1 and B NLSPs, but our results 
are easily extended to the general cases. 



2 Another common definition is Ynlsp = "nlsp/s, where s = (2ir 2 /45)g*sT 3 is the entropy density. In the 
era of NLSP decays to gravitinos, s ~ 7.0n^ G . 
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FIG. 1: NLSP lifetime in seconds (solid) and mass in GeV (dashed) in the (mg,<5m = wnlsp — 
m,Q — mz) plane for slepton and sneutrino NLSPs. 



A. Slepton NLSP 

The width for the decay of any sfermion to a gravitino is 

1 



r(/ - fG) 



48ttM2 m% 

* G 



1 

2 



(6) 



assuming the fermion mass is negligible. For Am = — itLq <C itLq, the sfermion decay 
lifetime is 



r(f -> fG) ~ 3.6 x 10 s s 



100 GeV 



Am 



m 



G 



1 TeV 



(7) 



The slepton lifetime and mass are given in the {tuq, 5m) plane in Fig. ^ where we have 
defined 

Sm = Am — mz = "Inlsp — m G ~ m z > (8) 

a useful measure of the kinematically available energy in three-body decays to be discussed 
below. 

For selectrons, the produced electron in these two-body decays immediately initiates an 
electromagnetic (EM) cascade, and so 



-°EM — 1 > e 



e 

EM 



2 2 

m| — m^ 



Imp 



(9) 



For smuons, the produced muon typically interacts with the background photons before 
decaying [l8|, and so 

2 2 

mf. — m>, , , 

$M = 2 • (10) 



-°EM — 



For staus, the resulting r sometimes decays into mesons, which could in principle induce 
hadronic cascades. As shown in Refs. however, for decay times r > 10 3 — 10 4 s, the 

hadronic interaction time of all pions and kaons is much longer than their decay time. The 
decays of staus therefore typically contribute only to EM cascades, and we assume this in 
the following analysis. In contrast to the selectron and smuon cases, however, on average, 
about half of the r energy is lost to neutrinos. We therefore have 

W ~1 e f ~ l m ?~ m 6 (11) 
em > 2 2m, • 1 j 

As noted in Sec. HJ three-body decays are also important when they are the leading 
contribution to hadronic cascades. They are therefore important for slepton NLSPs. The 
decays are those of Eq. (jSJ). The decay I — > IZG takes place through off-shell /, /, and \i 
and also through a four-point interaction. The three-body decay widths for sleptons have 
been discussed and presented in Ref. and we refer readers there for details. Given these 
decay widths, the hadronic branching fraction is 

r T(l^lZG)B^ d + r(f- vWG)BZ d + T{l->lqqG) 

had " r(i-ZG) ' { } 

for I = e, fx, f, where B^ ad ,B^ d « 0.7 are the Z and W hadronic branching fractions. 
r(Z — > uWG) = for purely right-handed sleptons. Below, we will consider cases in which 
the three-body decays are kinematically allowed. These decay modes may nevertheless 
become suppressed for Am ~ mz,mw- However, even for such small mass splittings, 
hadronic decays are still possible through higher order decays. With this in mind, we have 
included the four-body process T(l — ► IqqG). We have not calculated this width. However, 
we expect B(l — > IqqG) ~ 10~ 6 , and we take this value, which provides a lower limit on 

B l had is typically in the range — 10~ 5 , depending on the underlying scale and mass 
splitting. As the branching fraction may vary over a few orders of magnitude, variations in 
e l had are subdominant. We therefore take simply 

4ad = \ ( m l- m d) ( 13 ) 

in our analysis. 



B. Sneutrino NLSP 

The decay width and time for v — > vG are given in Eqs. © and (J7J), and plotted in 
Fig. ^ These two-body decays are essentially invisible and do not contribute to either EM 
or hadronic cascades. (We neglect the effects of neutrino thermalization through processes 
like t^eBG ~~ > ve.) The three-body decays are therefore even more important for sneutrinos 
than sleptons. These decays have also been discussed and presented in Ref. [l9j . For 
sneutrinos, we have 

r(z> -> vZG) + Y(y -> IWG) + T(u -> uffG) 
Bem ~ r(z> Z vG) (14) 
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e EM = ( 15 ) 



had Y{v^uG) { ' 

Chad = ^K- m c) , ( 17 ) 

for v = u e , Op, v T . The EM branching fraction is in fact slightly reduced by decays v — > z/ZG 
followed by Z — > i/p. We have neglected this effect. The four-body decay takes place through 
virtual neutralinos. We again assume B(p — > vqqG) ~ 10~ 6 , which provides a lower limit 

on #had- 

Note that, in contrast to the case of the slepton NLSP, the EM branching fraction is 
suppressed and of the same order as the hadronic branching fraction. In our analysis below 
we have included the EM constraint, but we find that it is so weak that it does not disfavor 
any of the parameter space appearing in figures below. The hadronic BBN constraint is so 
much stronger than the EM constraint at early times, however, that it is still important 
then, as we will see. 



C. Neutralino NLSP 



For neutralino NLSPs, the decay width to photons is 



v( _ \Nu cosOw + Nu sm6 w \ 2 m \ [ 



ml 



m 



x J 



3 r 



1 1 

1 + 3— f 

mr, 

A - 



(18) 



where \ = N n (-iB) + N 12 (-iW) + N 13 H d + N 14 H U . In the limit Am < m x , the decay 
lifetime is 



r(x -> 7^) ~ 2.3 x 10 7 s- 



cos 2 6 



w 



Nu cos6 w + N12 sin^l 2 



100 GeV 
Am 



(19) 



proportional to (Am) 3 and independent of the overall superpartner mass scale. 

This decay contributes only to EM energy. As noted in Sec. U the leading contribution 
to hadronic energy is from \ ~* ZG, hG. These decays produce EM energy for all possible 
Z and h decay modes (except Z — > vis), but they may also produce hadronic energy when 
followed by Z, h — > qq. The decay width to Z bosons is 3 



w ^ | — Nu sin 9 W + N12 cos 9w \ 2 m \ E - 



48ttM 2 

2 \ 2 / 
,2 



(20) 



Our expression for G in Eq. 1221) differs from the result of Ref. |3l[ in the sign of "12" in the second term. 
The authors of Ref. |3lJ used Eq. (4.31) of Ref. |32(, which contains a sign error. We have corrected for 
this error. We thank Y. Santoso and T. Moroi for helpful correspondence. 
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where 



F(m x ,m d ,m z ) 



1 - 



m 



tjIq + m z 



1/2 



m 



G(m x ,m 6 ,m z ) =3 + ^fl-12 + ^)+^f-^f(3-^f 



ms- ml 



m 



m 



mr, mr 



mr 



The decay width to the Higgs boson is 4 



r( x - hG) 



- N 13 sin a + JV14 cosal 2 m 



96ttM2 
2 / 



-F(m x ,m G ,m h ) 



G 



X. 



m 
m 



mi 



rn 



x 



where h = (—Hasina + if°cosa)/\/2, F is as given in Eq. (jUJ), and 

•J 

,2 



;;//, mj, mj, m 4 ^ mi 



m 



x 



rn 



x 



rn 



x 



X. A 



m 



M 3 + 2 I^ + 3 



m 



mf 



m. 



mf 



(21) 
(22) 



(23) 



(24) 



For the case of a Bino-like neutralino, the neutralino's mass and lifetime are given in the 



(ttIq, 5m) plane in Fig. El 



4 This result disagrees with the decay width given in Ref. |3l|. After cross-checking with the authors, they 
agree with our current results. We thank Y. Santoso and V. Spanos for helpful correspondence. 
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Given these two-body decay widths, the resulting values for the energy release parameters 

are 



Blu - 1 (25) 

^ = ^ <*> 
RX r( x -> ZG)Bl d + T(x - hG)Bt d + T( x -> ggG) 

m 2 x -ml + m\ h 
£had 2^ ' 

where B^ ad 0.9. For the three-body decay, we take T(x — > ggG) ~ 10~ 3 , which provides 
a lower bound on B£ ad when the two-body decays become kinematically suppressed. In 
addition, for m x — tuq < m,z, e£ad * s estimated to be §(m x — mg) in our analyses. 

We have neglected decays to the heavy Higgs bosons. When kinematically allowed, they 
will, of course, modify the branching fraction and energy release formulae above. The decay 
width to the heavy CP-even Higgs boson H is given by replacing rrih — > ran and —N13 sin a+ 
iVi4 cos a — > Ni 3 cos a + N u sin a in Eq. (J2*3|) . The decay width to the CP-odd Higgs boson 
A is given by replacing m h — > m^, — A/" 13 sina + iV 14 cosa — > N 13 sin j3 + A/" 14 cos/3 and, 
in the last line of Eq. (|23jh m x — > —m x , where the last transformation is required by the 
CP-odd nature of the A boson. 



III. THERMAL RELIC DENSITIES 



To determine the normalized NLSP number density Ynlsp of Eq. (jSJ) and also the resulting 
contribution of gravitinos to the current dark matter energy density, we assume that the 
NLSP freezes out with its thermal relic density. The super WIMP has no effect on the early 
thermal history of the Universe. The NLSP therefore freezes out as usual, with relic density 



given approximately by 33, 3 



th h 2 



1.1 x lO^GeV- 1 
y^Mpi cj(av) 



0.2 



15 



Xp 

30 



10 19 GeV 



Mi 



pi 



10~ 9 GeV" 2 
cj{av) 



(29) 



where xp = ?^nlsp/2f is the NLSP mass divided by the freeze out temperature Tp, g* is 
the effective number of massless degrees of freedom at freeze out, and (av) is the thermally- 
averaged NLSP annihilation cross section, and cj is 1 for S'-wave annihilation, 1/2 for P-wave 
annihilation. The energy release parameter Ynlsp is derived from this through 



NLSP 



^NLSPPc 
^NLSP^ G 



1.3 x 10 



-11 



TeV 

T^NLSP 



n 



th 

NLSP > 



and the gravitino relic density is given by 



n 6 h 2 



m 



"^NLSP 



"NLSP"- 



(30) 



(31) 



For the case of slepton NLSPs, the dominant annihilation channels are typically 11* — > 
77, jZ, ZZ through slepton exchange and 11 — > 11 through Bino exchange. For right-handed 
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sleptons, the thermally-averaged cross section near threshold may be approximated as [3j 



(<™)l 



Aira 2 

2 



+ 



cos 4 6Wm~ +m%) 2 



5.0 x 10" 10 C 



TeV 



GeV- 



(32) 



where m§, mj R are the Bino and slepton masses, respectively, and C is an 0(1) model- 
dependent constant. Here we have not included co- annihilation processes, which might be 
important if sleptons and, say, neutralinos are nearly mass degenerate. Using Eq. ({2Hj) and 
setting C — 1, the slepton thermal relic abundance is 



n\ h h 2 

lR 



0.2 



mr 



1 2 



TeV 



A similar analysis for the sneutrino NLSP case yields [2 



nfh 2 



0.06 



mo 



1 2 



TeV 



(33) 



(34) 



The thermal relic density of the sneutrino is typically smaller than that of right-handed 
sleptons because sneutrino annihilation is relatively efficient, taking place through weak 
SU(2) couplings, whereas the right-handed sleptons annihilate only through hypercharge 
couplings. 

For the neutralino NLSP case, the thermal relic density is very model-dependent. The 
annihilation cross section varies widely depending on the gaugino-Higgsino composition of 
the neutralino and the presence or absence of co-annihilation effects, and so depends on 
a large number of unknown supersymmetry parameters. Rather than constraining these 
parameters by working in a particular model, we adopt a simple scaling behavior based on 
some well-known results. In particular, we assume that the annihilation cross section scales 
as m~ 2 . To fix the constant of proportionality, we recall that in the "bulk" region of minimal 
supergravity, where the neutralino is Bino-like and there is no significant co-annihilation, 
the desired relic density is achieved for « 100 GeV. In the focus point (FP) region of 
minimal supergravity 3^, 37, 38[, where the neutralino is a Bino-Higgsino mixture [39], the 



neutralino mass may be much larger |40j, |41[ . If there are co-annihilation effects |42j, |4 



the neutralino mass may also be much higher To study the effect of having 

a heavier neutralino, we consider the mass m^ « 200 GeV as an example of these other 
possibilities. We therefore consider the range 



"bulk" 



A 



0.1 



m A 



1 2 



100 GeV 



to 



"focus point / co-annihilation" : 



Q^h 2 



0.1 



m A 



(35) 



(36) 



.200 GeVj 

Note that for similar NLSP masses, the thermal relic density is much higher in the 
neutralino case than in the slepton case. This is as expected, because the neutralino anni- 
hilation is dominantly P-wave because of the Majorana-ness of neutralinos, while slepton 
annihilation takes place in the S'-wave. In fact, for similar masses, one expects the slepton 
relic density to be suppressed relative to the neutralino relic density by a factor of roughly 
v 2 ~ 3/xp ~ 1/10. Given the approximations used, this is in reasonable quantitative 
agreement with the estimates of Eqs. (J3~3*|) - (}3l)j) . 
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IV. CONSTRAINTS 



A. Dark Matter Density 

An unambiguous and simple constraint on these scenarios is that the resulting gravitino 
energy density should not be greater than the observed dark matter density. This constraint 
may be avoided if there is significant entropy production between NLSP freeze out and now. 
However, assuming such new physics is counter to our goal of evaluating the gravitino LSP 
possibility in the simplest possible cosmology, and so we require 

n~ h 2 = -!^n$ LSP h 2 < 0.11 . (37) 

m NLSP 

As evident from Eqs. (JTj) and (fT9")) and Figs. [T] and |2l the typical decay times are t ~ 10 4 — 
10 8 s. This is the natural decay time of a particle with weak-scale mass that decays through 
gravitational interactions. There are therefore additional constraints, most importantly from 
bounds on EM energy release from cosmic microwave background (CMB) \i distortions and 
from bounds on both EM and hadronic energy release from BBN light element abundances. 



B. Cosmic Microwave Background 

The CMB constraint is fairly straightforward to understand. The CMB photon energy 
distribution at times t ^5 10 8 s may be parameterized as 

f ^ = e E/(kT)+n _ I ' ( 38 ) 

with chemical potential /i. For early decays, EM cascades are completely thermalized 
through energy-changing processes je~ —>■ and number-changing interactions, such 
as eX — > eX'-f, where X is an ion, and double Compton scattering 7e~ — > 77c". The 
resulting distribution is therefore Planckian, with \x = 0. For decay times in the window 
of interest, however, the number-changing processes may be inefficient. In this case, the 
spectrum cannot relax to a distribution determined by only one parameter, the tempera- 
ture T. It therefore relaxes to statistical but not thermodynamic equilibrium, resulting in a 
Bose-Einstein distribution function with \x 7^ 0. 

The value of the chemical potential /1 may be approximated for small energy releases by 
analytic expressions given in Ref. |47| . These have been updated with current cosmological 
parameters in Ref. . We will apply the current constraint jH, Hi| 

< 9 x 1(T 5 . (39) 



C. Big Bang Nucleosynthesis 



The BBN constraints are more complicated and more ambiguous. Constraints on EM 

51, 53[. Most recently, EM constraints 
(but not hadronic constraints) have been considered in Ref. j^ij and these were used in the 
previous analyses of Refs. [13, EH- Here we include contours corresponding to the most 



energy release have been studied in 0, H, IHol 

i ts) hav< 

stringent constraint from that analysis, the deuterium bound 

1.3 x 1(T 5 < D/H < 5.3x10" 



(40) 
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2.4 x 1(T 5 


< D/H 


< 


3.2 x 10 




3 He/D 


< 


1.13 




6 Li/H 


< 


6.1 x 10 


0.228 


< Y p 


< 


0.248 . 



to facilitate comparison with previous results. 

More recently, both EM and hadronic energy releases have been bounded in the analysis 
of Ref. [23j]. Of the constraints imposed there, the most relevant for us are the 2a bounds 

(41) 
(42) 

11 (43) 

(44) 

The statistical and systematic errors have been combined in quadrature for the 4 He (Y p ) 
constraint, and the 6 Li/H result is obtained by combining the 95% confidence level (CL) 
constraints 6 Li/ 7 Li = 0.05 ± 0.02 and 7 Li/H = 2.2±f;| x 10" 10 Q. 

As is evident, the later analysis is much less conservative. First, it assumes a significantly 
more stringent D bound. Measurements of primordial D have long been considered by many 
to be the most reliable baryometers. There is also now impressive concordance between the 
baryon number determinations from D and CMB measurements, which further supports the 
narrow range of D/H given in Eq. (|41j) . At the same time, existing discrepancies between 
standard BBN predictions and observations in other elements may indicate that caution is 
still needed in interpreting the D bound. In particular, if these discrepancies are indications 
of new physics, the required new physics is also likely to distort the D abundance, since the 
D binding energy is so small. 

The analysis of Ref. [23| also includes stringent constraints from 3 He/D and 6 Li/H. If 
taken at face value, these additional bounds in fact provide some of the most stringent 
bounds on the gravitino LSP scenario. The 3 He/D bound is the strongest constraint on EM 
energy and provides the strongest constraint for NLSP decay times r > 10 7 s, while the 6 Li/H 
constraint on hadronic energy release provides the strongest constraint on earlier decay times. 
At the same time, it is important to bear in mind that these constraints are on less sure 
footing than the D constraints. For 3 He, 3 He/H suffers from uncertainties in chemical/stellar 
evolutio n |56j| . Although 3 He/D has been proposed as an alternative tool to constrain new 
physics |5lL l57l|. present evaluations exist only in the Sun [HH, and the determination of 
primordial 3 He/D requires a rather involved extrapolation of these results. 6 Li has also 
been proposed as a promising probe of new physics. However, after WMAP, there is a clear 
discrepancy between standard BBN predictions and the observations of 7 Li 
consistency possible only if systematic uncertainties have been underestimated |62 
calls the status of 6 Li into question, as direct observations of 6 Li/H are difficult, and so the 
upper limit on 6 Li/H is usually derived from bounds on 6 Li/ 7 Li. In fact, the current status 
of 6 Li and 7 Li may also be taken as evidence for new particle physics 6^, 64j . 




In light of all of these comments, we present constraint contours from all of these data, 
including both D constraints, to show the (strong) effect of varying BBN assumptions. We 
consider regions of parameter space that violate the conservative constraint of Eq. ()40|) to be 
excluded, but we consider regions that violate only Eqs. (j4*Tj) - (j4^j) to be at most disfavored, 
but not necessarily excluded, given the significant ambiguities noted above. 

There are subtleties in importing the constraints of Eqs. (|4ip -(|43 |) to the present analysis. 
When including both EM and hadronic energy release, there is the possibility of cancella- 
tions. In addition, although the EM constraint depends essentially only on £ EM of Eq. 
the hadronic constraint may depend, in principle, ehad and -Bhad^NLSP separately, and results 
are presented only for a few values of ehad- In practice, however, the cancellations occur 
only in rather special cases for particular energy release time. In addition, in supergravity 
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with a gravitino LSP, the NLSP lifetime is usually larger than 150 sec, and so the hadronic 
constraint depends to a good approximation on £had = Chad -Shad ^nlsp only. 5 We therefore 
impose constraints on £had and impose the constraints on EM and hadronic energy release 
separately, ignoring the possibility of cancellations. 

V. RESULTS 

We have now determined the energy release parameters B\ and e\ in Sec. HU and Inlsp 
in Sec. 11111 We may now compare these to the constraints of Sec. I1VI to determine what 
combinations of gravitino mass and NLSP mass are excluded, disfavored, and allowed for 
various NLSP possibilities. We will present results in the {m^, 5m) plane, where 5m = 
^nlsp ~~ m G ~ m z- We consider only 5m > 0, so three-body decays are therefore always 
kinematically possible. Of course, they are highly suppressed for small 5m. 

A. Slepton NLSP 

We begin with the stau NLSP scenario. We assume the stau is right-handed. Neutralino 
and chargino parameters enter in the three-body decay widths. We take fi = M 2 = 2Mi = 
4mf H and tan f3 = 10. 

The results are presented in Fig. EJ To understand these results, it may be helpful to refer 
to the mass and lifetime contours of Fig.^ Note that, given the definition of 5m, mr R > m z 
in the entire plane. The current limit on a metastable stau from LEP is mf R > 99 GeV [65], 
and so excludes a small portion of the lower lefthand corner. 

The shaded regions are excluded. Given the scaling ^nlsp ^ m NLSP f° r the thermal 
relic density, the dark matter density implies an upper bound on the product m^spm^,, 
excluding f and gravitino masses ~ 1 TeV. The constraint is relatively mild, because staus 
annihilate efficiently through S'-wave processes. The other shaded region is excluded by the 
absence of CMB /j, distortions. This provides a more stringent constraint than Qq for decay 
times t > 10 7 s, when the decay products are produced too late to be thermalized. 

The BBN sensitivity contours divided into those from D and 4 He, which are probably 
the most reliable (left panel), and those from 3 He and 6 Li, which are on less sure footing, 
given the discussion of Sec. II V CI For the D and 4 He results, we present results given 
the conservative bound of Eq. (|40j) on EM cascades (EMI), the more aggressive bound of 
Eqs. flUJ) and flSJ) on EM energy (EM2), and the bound of Eqs. (JHJ) and flEJ) on hadronic 
cascades (had). The EMI contour lies completely in the CMB-excluded region. Although 
BBN constraints are often assumed to be the leading constraint on late decays, we see that 
the CMB spectrum is now know to be Planckian to such high precision that the CMB 
constraint is competitive with the leading BBN constraints. At the same time, we see that 
the strength of the EM constraints in constraining gravitino LSP parameter space depends 
sensitively on how one interprets the BBN data. Adopting the more stringent EM2 contour, 
we find that bounds on are improved by about an order of magnitude. 

Our analysis includes hadronic bounds on the gravitino LSP scenario for the first time. 
From Fig. El we see that the hadronic constraint is the leading constraint for relatively 



5 The constraint on £had varies only by a factor of 2 for ehad between 100 GeV and 1 TeV [5 



13 



1000 
500 

100 b- 



% 50 

a 



10 
5 



D / H+Y P had ^ 



EM2, 



CMB 




1000 
500 

100 b- 



% 50 



10 
5 



3 He/D EM 




5 10 50 100 

mg (GeV) 



5001000 



5 10 50 100 5001000 

mg (GeV) 



FIG. 3: Excluded and allowed regions of the (mg, 8m = mNLSP — m Q — m z) parameter space in the 
gravitino LSP scenario, assuming a tr NLSP that freezes out with thermal relic density given by 
Eq. (|33|). The light (yellow) shaded region is excluded by the overclosure constraint fl^h? < 0.11, 
and the medium (green) shaded region is excluded by the absence of CMB \i distortions. BBN 
is sensitive to the regions to the right of the labeled contours. Left: Regions probed by D and 
4 He, assuming the conservative result of Eq. (|40jl (EMI), and the more stringent constraints of 
Eqs. Q41|) and (|44|1 (EM2 and had). The dotted line denotes the region where cancellation between 
D destruction and creation via late time EM injection is possible 27], while 7 Li is reduced to 
the observed value by the late NLSP decays [54j]. Right: Regions probed by 3 He/D (EM), 6 Li/H 
(had) 0, and 6 Li/H (EM) [HQ. 



early decays. Recall that sleptons produce hadronic energy only in three-body decays, and 
so the hadronic energy release is suppressed by factors of ~ 10~ 3 relative to EM energy. 
Nevertheless, hadronic decay products are so lethal to light elements that the hadronic 
constraints are the most stringent constraint in parts of parameter space. Note also that the 
part of parameter space in which hadronic constraints are most important is where and 
ttlq are both in the hundreds of GeV, the most natural region for weak-scale supergravity. 
We conclude that hadronic constraints and three-body decays must be taken into account 
to establish the viability of any gravitino LSP scenario. 

In the right panel of Fig. EJ we include the sensitivity contours of 3 He and 6 Li. We see 
that, taken literally, the constraint on 3 He provides the most stringent constraint on late 
decays (through its limits on EM energy) and 6 Li provides the leading constraint on early 
decays (through its limits on hadronic energy). Of course, given the ambiguities discussed 
in Sec. II V CI we do not consider these contours to be exclusion contours. These sensitivity 
contours are of interest, however, as, if colliders measure the superpartner parameters to be in 
the regions to the right of these contours (as we will discuss in Sec. IV1J). these measurements 
will have important implications for BBN. 

Last, we discuss the dark matter implications of these results. At the boundary of the 
region excluded by Qq, the light (yellow) shaded region, gravitino superWIMPs account for 
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FIG. 4: As in Fig. |HJ but assuming a sneutrino NLSP that freezes out with thermal relic density 
given by Eq. (|33)). 

all of dark matter. Taking the EMI and hadronic constraints from D and 4 He, we see that 
this possibility is indeed viable. If, however, the various BBN anomalies are resolved and 
the EM2 and 3 He and 6 Li contours may be considered as exclusion contours, the possibility 
that superWIMP gravitinos form all of dark matter may be excluded. However, given the 
current status of BBN, we find such conclusions premature. 

Finally, we should remind the reader that we have assumed a particular thermal relic 
density and particular neutralino mass parameters, which enter the three-body branching 
ratios. All of the contours above will shift if there are significant deviations in these assump- 
tions. A complete analysis of all of these variations is, however, beyond the scope of this 
work. 

B. Sneutrino NLSP 

As discussed in Sec. I11B1 if the NLSP is a sneutrino, two-body decays are essentially 
invisible, and so inclusion of three-body decays is essential to determine the viable parameter 
space. Because the hadronic constraints are some much stronger than the EM constraints, we 
may focus on them only. Again, neutralino and chargino parameters enter in the three-body 
decay widths, and we assume \x = Mi = 2M\ = 4m„ and tan/3 = 10. 

The results are presented in Fig. 0] Sneutrinos annihilate through S"-wave processes even 
more efficiently than sleptons, as can be seen by comparing Eqs. ()33|) and (|34|). The dark 
matter density bound is therefore weaker, and is in fact pushed to the right; it does not 
appear in the plotted plane. The CMB constraint, previously so stringent, is also absent, of 
course, as it constrains EM energy, and the BBN constraints on EM energy are also absent. 

The remaining constraints are therefore only the hadronic BBN constraints. These are 
stringent for early decays, that is, large 5m. The more reliable D and 4 He constraints disfavor 
5m > 300 GeV, while 6 Li (had) is sensitive to 5m > 200 GeV. It is rather remarkable 



15 




FIG. 5: As in Fig. OJ but assuming a Bino NLSP that freezes out with the "bulk" thermal relic 
density given by Eq. (|35]> . 



that the sneutrino NLSP case is so tightly constrained, given its invisible dominant decay 
mode. At the same time, the scenario is perfectly viable for natural weak-scale supergravity 
parameters. Note that gravitino superWIMP dark matter is also viable for ~ 1 TeV 
and 8m < 300 GeV. 



C. Bino NLSP 

Finally, we turn to the case of the Bino NLSP. The results are presented in Fig. El for the 
case where the Bino thermal relic density is as in the bulk region of minimal supergravity 
(Eq. (jHSJ); and in Fig. [H] for the case where the Bino relic density is degraded by a factor of 
4 (Eq. (jUS))), as might be the case if there there are additional effects as may be found in 
the focus point or co-annihilation regions of minimal supergravity. Note that the tUq lower 
limit has been extended to much lower masses than in the slepton and sneutrino figures. 
The mass and lifetime contours of Fig. Elmay be helpful in understanding these results. 

The CMB and EM BBN bounds are roughly similar to those in the slepton NLSP case. 
However, for the other bounds, there are important changes. Because neutralino annihila- 
tion is P-wave suppressed, the dark matter density limit is much more stringent, excluding 
gravitino masses above 100 GeV and 200 GeV in the "bulk" and "FP/coann" cases, respec- 
tively. 

Even more striking, the hadronic BBN bounds become much more stringent. This is 
expected — Bino NLSP decays contribute to hadronic energy at the two-body level through 
decays B — > ZG. The branching fractions for EM and hadronic decays are therefore not too 
different, and the extreme stringency of the hadronic constraints makes them the dominant 
bound. Taking only the relatively reliable D and 4 He bounds, we find that decay times 
t > 10 3 s are disfavored, excluding almost all gravitino masses tjiq ^ 100 MeV. In the 
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FIG. 6: As in Fig. G3 but assuming a Bino NLSP that freezes out with the "focus point/co- 
annihilation" thermal relic density given by Eq. (|36|) . In the left panel, the region between the 
"had" lines is disfavored. 



"FP/coann" case, for 5m < O(10 GeV), the hadronic constraints exclude 100 MeV < w,q ^ 
500 MeV, but tuq ~ 1 GeV is again allowed. In the allowed region, the decay times are so 
long t > 10 5 s that the hadronic constraints become less stringent. However, for such long 
decays, the EM constraints are stringent. As can be seen in the figures, the EM constraints 
disfavor or exclude this island of parameter space, depending on how conservative one's 
interpretation of the EM BBN constraints is. 

We therefore conclude that hadronic BBN constraints essentially exclude supergravity 
with a gravitino LSP and a Bino NSLP when the decay channel B — > ZG is open. The 
Bino NLSP scenario is viable only if the Bino and gravitino masses are degenerate enough to 
suppress this decay mode, if tuq is below 10 MeV, a rather unnatural value of conventional 
supergravity, or, possibly, if there are extremely co-annihilation effects which suppress the 
thermal relic density even more than in our "FP/coann" example. If the neutralino NLSP is 
not pure Bino, there are additional possibilities. For example, as noted in Ref. 3], photino 



NLSPs may be viable, as they contribute to hadronic cascades only through three-body 
decays. Of course, from the high energy viewpoint, a photino-like neutralino is unmotivated. 
More likely is the case of a Higgsino-gaugino neutralino, for which the thermal relic density 
may also be greatly suppressed, but the hadronic constraints are stronger, since T(x — > hG) 
is larger. A detailed examination of such focus point or co-annihilation cases would be 
interesting, but is beyond the scope of this study. 



VI. IMPLICATIONS FOR COLLIDER PHYSICS 



The possibility of a gravitino LSP in supergravity has rich implications for current and 
future colliders. These implications depend crucially on whether the NLSP is a slepton, 
sneutrino, or neutralino. In all cases, however, given NLSP rest lifetimes of 10 4 — 10 8 s, the 
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FIG. 7: Excluded and allowed regions for the gravitino LSP scenario, assuming a tr NLSP, as 
in the left panel of Fig. |3J but now in the (m^,mf R ) plane. In the dark (red) shaded region the 
gravitino is not the LSP. All other shaded regions and contours are as in Fig. |21 

typical NLSP decay lengths are enormous relative to collider detectors, and so these NLSPs 
are essentially stable as far as colliders are concerned. 

A. Slepton NLSP 

We begin by discussing the slepton NLSP case. This possibility is very natural from the 
point of view of high-energy frameworks. Given simple boundary conditions at the grand 
unified scale, for example, and evolving these to the weak scale, right-handed sleptons, in 
particular, right-handed staus, often emerge as the lightest standard model superpartner. 
In conventional studies of these high energy frameworks, such regions of parameter space 
are excluded by bounds from searches for charged massive stable particles in sea water. In 
the gravitino LSP scenario, however, the lightest slepton is metastable but not absolutely 
stable, and so these bounds do not apply. 

The excluded and allowed regions of parameter space are presented again in Fig. [7[ but 
now in the (m^, Trif a ) plane, which is more convenient for inferring implications for colliders. 
We see that, allowing the gravitino to be as light as 10 GeV, all weak-scale stau masses are 
allowed. Staus may therefore be within reach of the LHC and even of the first stage of a 
linear collider. For heavier gravitino masses, the allowed stau mass range becomes more 
narrow. Neglecting the aggressive EM2 bound, we see that all dark matter may be in the 
form of gravitino superWIMPs if staus have masses m^ R > 1 TeV. In this case, direct 
stau production is beyond the range of the LHC and linear collider, but staus may still be 
produced in the cascade decays of squarks and gluinos. 

At hadron colliders, sleptons can be pair-produced through the Drell-Yan processes 

qq' — > W* — > l L H (45) 
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qq -> Z*,7* -> Z L Z L , [r/r, z/ l z/ l . (46) 

The cross sections for such processes are determined by the slepton masses, with very little 
other model dependence. 6 These Drell-Yan cross sections have been studied in detail [H^, 



671 . l68j . including the leading QCD corrections. For the Tevatron with y/s = 2 TeV and 
raj = 100 GeV, the cross sections for IrIr, IlIl (v L i) L ), and XiPl are about 10 fb, 30 
fb, and 100 fb, respectively. These cross sections drop quickly for heavier sleptons. For 
mi = 200 GeV, they are reduced by more than an order of magnitude, and sleptons with 
mass above around 250 GeV will be beyond the reach of the Tevatron. 

For the LHC with <Js = 14 TeV, the Drell-Yan cross sections are about 10 times bigger. 



Sleptons may also be produced via weak boson fusion 

qq' -> qq'VV -> qq'tl . (47) 

This cross section decreases much more slowly with increasing slepton mass than the Drell- 
Yan cross section, and the weak boson fusion cross section dominates for mj > 200—300 GeV. 
At the LHC, hundreds to thousands of sleptons could be produced. 

Metastable sleptons will appear as charged tracks in the tracking chamber with little 
calorimeter activity. Eventually they will hit the muon chambers and so look muon-like. 
However, given their large mass, such sleptons may be non-relativistic. They can therefore 
be highly ionizing, allowing one to distinguish them from genuine muons. In addition, time- 
of-flight information could be used to detect a slow moving particle. Suc h sig nals are almost 
background free, providing the potential for a spectacular signature 

Searches for metastable sleptons have been motivated previously by the existence of such 
particles in gauge- mediated supersymmetry breaking models 6f| E3, El, Eil EEf • No signals 
have been found at Tevatron Runs I and II 76[ and LEP (65 1. The most stringent current 



bound is mj > 99 GeV at 95% CL from LEP searches at center-of-mass energies up to 208 
GeV. The prospects for a full Tevatron Run II have been investigated in Ref. |66j], where 
appropriate cuts in slepton velocity and pseudorapidity rj have been included to eliminate 
the background. Requiring 5 or more events for a signal, the estimated reach in right handed 
slepton mass is about 110 GeV, 180 GeV, and 230 GeV for integrated luminosities of 2, 10, 
and 30 fb -1 , respectively. The discovery reach of the LHC has also been considered [T7| . 
For one year at the design luminosity of 100 fb" 1 , metastable sleptons with mass up to 700 
GeV could be discovered. 

As noted above, metastable sleptons are possible in both the high-scale supersymmetry 
breaking scenarios discussed here, and in low-scale supersymmetry breaking models, such 
as those with gauge-mediated supersymmetry breaking. In the gauge-mediated scenarios, 
however, the gravitino mass is much lighter, around the keV scale. It is possible that these 
cases may be distinguished cosmologically. Alternatively, direct collider searches for other 
supersymmetric particles and the measurement of their mass spectra will provide additional 
means for distinguishing these possibilities. Finally, the slepton lifetimes in gauge-mediated 
models, although long on collider detector scales, are much shorter than in the superWIMP 
scenarios discussed here, and this may be distinguished experimentally, providing an unam- 
biguous determination |78j |. 



6 Much larger cross sections may result from sleptons produced in cascade decays of gluinos and squarks, 
but the details of these processes are highly model-dependent. 
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FIG. 8: Excluded and allowed regions for the gravitino LSP scenario, assuming a v NLSP, as in 
Fig. 0J but now in the (m^mp) plane. In the dark (red) shaded region the gravitino is not the 
LSP. The D + 4 He contour is as in Fig. 0] 

B. Sneutrino and Neutralino NLSP 

The cases of a gravitino LSP with either a sneutrino or neutralino NLSP are qualitatively 
different from the slepton NLSP case. The allowed regions of parameter space are given in 
Figs. |H1 and El In both cases, the metastable NLSP will pass through detectors, resulting in 
missing energy signatures topologically identical to the conventional missing energy signal 
of supersymmetry. There are four cases to distinguish: the lightest standard model super- 
partner may be either a sneutrino or a neutralino, and this particle may either decay to a 
gravitino or not. 

The sneutrino and neutralino cases may be distinguished by precision supersymme- 
try studies at colliders. For example, in e + e~ collisions, the signatures of slepton pair 
production in the sneutrino scenario (for example, e + e~ — > 11 — > luuqq'u) are iden- 
tical to the signatures of charging production in the neutralino scenario (for example, 
e + e~ — > x + X~ ~~ > lvx°qq'x°) [79|- However, these possibilities may be distinguished easily 
through angular distributions at a linear collider, and possibly also at the LHC. 

Determining whether the sneutrino or neutralino eventually decays to a gravitino may be 
more difficult. In the sneutrino case, the working assumption would be that the gravitino 
is the LSP — the sneutrino itself is disfavored as a dark matter candidate, because, in 
the natural region of parameter space, it predicts dark matter signals that have not been 
seen. On the other hand, in the neutralino case, the working assumption would be that 
the neutralino is stable — as we have seen, decays to the gravitino are highly constrained 
by hadronic BBN bounds. Further information will be provided by high sensitivity collider 
and astrophysical experiments. For example, a positive detection in dark matter search 
experiments would eliminate the possibility of a gravitino LSP. On the other hand, the 
precise determination of supersymmetry parameters will make possible the determination of 
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FIG. 9: Excluded and allowed regions for the gravitino LSP scenario, assuming a B NLSP in the 
(niQ,m^) plane. In the (dark) red shaded region the gravitino is not the LSP. All other shaded 
regions and contours are as in the left panel of Fig. 03 (left) and as in the left panel of Fig. |H1 (right) . 

the thermal relic abundance of the lightest standard model superpartner. This will favor the 
gravitino LSP scenario if this thermal relic density is larger than the observed dark matter 
density. 

The current experimental limit on the sneutrino mass, assuming three degenerate families, 
is nit > 44.6 GeV at 95% CL from limits on the invisible decay width of the Z |8(|. The 
metastable sneutrino scenario has not been well-studied. The collider signals and detector 
reaches depend crucially on the identity of the second lightest sparticle (which decays into the 
sneutrino), the accompanying decay products, experimental cuts and details of the detector. 
A detailed study of the collider phenomenology of the sneutrino scenario will appear in a 
future work. 

VII. SUMMARY AND OUTLOOK 

In this paper, we have determined the viability of supergravity scenarios in which the 
gravitino is the LSP. We have considered the three possibilities in which the NLSP is a 
slepton, a sneutrino, and a neutralino. In each case, we determined the branching fractions 
of the leading two- and three-body decays, and applied constraints from the dark matter 
density, CMB, and both EM and hadronic BBN bounds. We found that the hadronic BBN 
constraints, previously neglected, are extremely important, providing the most stringent 
limits in natural regions of parameter space. 

The gravitino LSP scenario opens up many connections between particle physics and 
cosmology. Consider, for example, the slepton NLSP scenario. At colliders, it may be 
possible to collect the less energetic metastable sleptons in a detector and monitor this 
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detector for slepton decays. By measuring the decay time distribution and the energy of 
each produced lepton, one could independently determine both the gravitino mass and the 
reduced Planck mass M* 2^, 3^ • 



A measurement of the gravitino mass determines the scale of supersymmetry breaking, 
with implications for model building and dark energy. At the same time, such a measurement 
would determine a particular place in the (rriQ, wnlsp) plane. Colliders would therefore 
shed light on the possible role of such new physics in BBN, and more generally in the 
thermal history of the Universe after NLSP freeze out at temperatures of around 10 GeV. 
For example, the fi parameter sensitivity of Eq. (l3"U{) may be improved by two orders of 
magnitude in the future by the DIMES mission j8l|. Such improvement would extend the 
CMB sensitivity contour significantly. The measurement of a fi distortion consistent with 
the determination of (m^, m-NLSp) would provide a striking confirmation of the underlying 
gravitino scenario. 

The measurement of M* would provide a precision test of the supersymmetry predictions 
relating the properties of gravitinos to those of gravitons, and also provide the first direct 
measurement of the Planck scale on microscopic scales jiil. l3(il] . The crucial question is the 
feasibility of collecting a sizable sample of metastable NLSPs. There has been an earlier 
study on the collection of very long lived heavy charged leptons [7l|. A more detailed 
analysis of the trapping of sleptons at future colliders is also now under study |7^ |. 
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